Heavy quarkonium production: 
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Two novel QCD effects, double color filtering and mutual boosting of the saturation scales in colliding 
nuclei, affect the transparency of the nuclei for quark dipoles in comparison with proton-nucleus 
collisions. The former effect increases the survival probability of the dipoles, since color filtering in 
one nucleus makes the other one more transparent. The second effect acts in the opposite direction 
and is stronger, it makes the colliding nuclei more opaque than in the case of pA collisions. As a 
result of parton saturation in nuclei the effective scale is shifted upwards, what leads to an increase 
of the gluon density at small x. This in turn leads to a stronger transverse momentum broadening 
in AA compared with pA collisions, i.e. to an additional growth of the saturation momentum. Such 
a mutual boosting leads to a system of reciprocity equations, which result in a saturation scale, a 
few times higher in AA than in pA collisions at the energies of LHC. Since the dipole cross section is 
proportional to the saturation momentum squared, the nuclei become much more opaque for dipoles 
in AA than in pA collisions. For the same reason gluon shadowing turns out to be boosted to a 
larger magnitude compared with the product of the gluon shadowing factors in each of the colliding 
nuclei. All these effects make it more difficult to establish a baseline for anomalous J/$ suppression 
in heavy ion collisions at high energies. 

PACS numbers: 24.85.+p, 25.75.Bh, 25.75.Cj, 14.40.Pq 



I. INTRODUCTION 

Nuclear suppression of heavy quarkonia is usually con- 
sidered as a sensitive hard probe for the properties of the 
short-living medium produced in heavy ion collisions [J- 
3] . The main challenge is to discriminate between initial 
state interactions (ISI) , usually identified as cold nuclear 
matter effects, and final state interaction (FSI), which is 
related to attenuation of the produced quarkonium in the 
dense matter created in the nuclear collision. While the 
latter is the main goal of the study, experimental infor- 
mation fully depends on how well we understand the ISI 
contribution. 

Naturally, the ISI dynamics should be studied in 
proton-nucleus collision, where no dense matter is ex- 
pected to be produced. The next step, the extrapolation 
of the results to nuclear collisions, is not that easy. Usu- 
ally it is done in an oversimplified manner, assuming that 
a cc pair is produced momentarily inside the nucleus and 
then attenuates with an unknown absorption cross sec- 
tion on the way out of the nucleus. The cc break-up cross 
section is fitted to pA data and used to predict the ISI 
effects in AA collisions. As is shown below, the break-up 
cross section is not constant, but steeply rises with the cc 
energy, besides, it is well known from HERA data. The 
most appealing oversimplification, especially at the high 
energies of RHIC and LHC, is the instantaneous produc- 
tion of cc. In reality the production time ranges from tens 
to thousands fermi at these energies. The long produc- 
tion time leads to shadowing which is a high twist effect 



suppressed by the quark mass squared. Nevertheless, we 
find that this shadowing effect is stronger than the lead- 
ing twist gluon shadowing, which produces a rather mild 
suppression even at LHC in the central rapidity region. 

Even if the shadowing and break-up effects in pA colli- 
sions are under control , still one cannot predict the cold 
nuclear effects in AA collisions in a model independent 
way. The two major phenomena considered in this paper 
make this impossible. The first one is double color filter- 
ing of dipoles propagating simultaneously through two 
nuclei. As is demonstrated below, the survival probabil- 
ity of a dipole in AA collisions is higher than the product 
of that in each of the colliding nuclei. 

Another effect changing the properties of the cold nu- 
clear medium in A A collisions, is the boosting of the satu- 
ration scales in the nuclei due to mutual multiple interac- 
tions in the nuclei. The effect is especially strong at the 
energies of LHC. The nuclear medium becomes several 
times more opaque for dipoles compared with its trans- 
parency in pA collisions. These two effects essentially 
modify the effects of the ISI stage of heavy ion collisions, 
which is usually considered as the baseline for search for 
" anomalous" nuclear suppression of heavy quarkonia. 

All calculations throughout this paper are performed 
in so called "frozen" approximation, neglecting the size 
fluctuations of dipoles during their propagation through 
nuclei. This approximation is accurate provided that the 
coherence length for heavy quark production substan- 
tially exceeds the size of the nuclei. This condition is 
satisfied at the energies of RHIC and above. 



II. PROTON-NUCLEUS COLLISIONS 

First of all one should make sure that the nuclear ef- 
fects for heavy quarkonium production in pA collisions 
are understood, as well as the absolute magnitude of the 
cross section. The cross sections of \ and J/\E' produc- 
tion in pp collisions was successfully reproduced in the 
color singlet model in refs. @ and @ respectively. What 
has been missed in most of current analyses of J/\& pro- 
duction in pA collisions, are the coherence effects related 
to the long time scale of charm production, associated 
with propagation of a color-octet cc pair through the nu- 
cleus prior the production of a colorless dipole [5, 7] . This 
stage results in the heavy quark shadowing, which is of 
the same order as the effect of break-up of the colorless 
dipole. Although both are the high twist effects, quanti- 
tatively they are more important for currently available 
data, than the suppression caused by the leading twist 
gluon shadowing. The latter is frequently miscalculated 
basing on some of nuclear gluon PDFs, which rely on ei- 
ther ad hoc or incorrect assumptions (see discussion in 



A. High twist heavy quark shadowing 

Although the proper time of charm production is short, 
t* ~ l/(2m c ), this time rises with J/ty energy linearly, 
in the rest frame of the nucleus, 



2E 



(1) 



Thus, if the energy of the produced J/^S is sufficiently 
high, E > 25( GeV) x L( fm), effects of coherence become 
significant. At the energy of RHIC, y/s = 200 GeV, and 
positive rapidities, t c > 12 fm. For T production t c at 
RHIC is rather short, but becomes much longer than the 
nuclear size at the energies of LHC. In what follows we 
assume the coherence time to be much longer than the 
nuclear size, unless otherwise specified. 

The nuclear suppression caused by coherence can be 
interpreted as high twist shadowing in the process of cc 
pair production by a projectile gluon. The cc is produced 
coherently in multiple interactions of the projectile gluon 
with target nucleons. If the J/*f? energy is sufficiently 
high, so that t c 3> R A , one can neglect the dipole size 
fluctuations during propagation through the nucleus, i.e. 
treat the dipoles being " frozen" by Lorentz time dilation. 

Since the production amplitude is convoluted with the 
charmonium wave function, one can assume with good 
accuracy an equal sharing of the total longitudinal mo- 
mentum between c and c. In what follows we rely on the 
saturated parametrization of the dipole cross section Q 



Vcc(r T ,x 2 )=a [l - e-^A-o 2 ^)" , 
where <jq = 23.03mb; rg(x2) = 0.4fm x (a^/xo) 



(2) 



x = 3.04 x 10- 4 ; x 2 = e~v ^(M| c ) + (p 2 T ) j yfl. The cc 

invariant mass distribution predicted by the color singlet 
model leads to (Mf c ) = 2Mj /iS 0. The measured (p^) = 

4 GeV 2 . 

The amplitude of cc production at a point with impact 
parameter b and longitudinal coordinate z inside the nu- 
cleus, averaged over the dipole size, reads [5J, 



S pA (b,z) = / d 2 r T W- cc (r T ) 



(3) 



x exp 



,(r T )T-(b,z) - -a ec (r T )T + (b, z) 



Here T_(b,z) = f'^dz' p A (b, z'); T+(b,z) = T A (b) - 
T_(6,z), and T A (b) = T_(&,oo). According to [| 
shadowing for cc production over the nuclear thickness 
T_(6, z) occurs with the shadowing cross section cor- 
responding to a 3-body dipole, gluon and cc, which 
for equal momenta of c and c equals to ascgirr) = 
Z a cc(r T /2) - |cr ec (r T ). 

The survival probability amplitude Eq. (J3J should be 
squared and integrated over the coordinate of the pro- 
duction point. Then the nuclear ratio reads, 



dz \S pA (b,z)\ x 



(4) 



We rely on the form of nuclear effects for production of 
the P-wave charmonium \2 derived in The mecha- 
nism of J/^ production is more complicated, but the gen- 
eral structure of shadowing corrections should be similar, 
because the overlap of the initial narrow size distribution 
of the perturbatively produced cc pair with the large size 
wave function of the final charmonium is dominated by 
the small initial size ~ l/m c , whether the charmonium 
is x ° r J/^- In addition, \ decays contribute about 30% 
to J/ 1 ^ production. 

In the regime of long t c ^> R A the weight factor in ([3]) 
has the form [a], 



W 5c (r T ) oc K (m c r T ) r 2 T ^j/*(r T ), 



(5) 



where one factor rx comes from the amplitude of cc 
production, and another one either from the amplitude 
of gluon radiation in the case J/\& production, or from 
the radial wave function of X2- Further, we assume the 
wave function to have the oscillatory form, ^ j/^{tt) oc 
cxp[-r^/2(r 2 J/qi )} with (r 2 ^) = 2/m c u, m c = 1.5 GeV, 
and lu = 300 MeV 0. 

To simplify the calculations we use for ^LfVx) = 
r 2 ^ Ko(rTm c ) the parametrization proposed in [8|, 



*i„(rr) w const x [e^/' 1 " - e" r ' /b 



(6) 



0.144. 



where a c — 0.496 fm, b c = 0.11 fm. The comparison 
of the two functions Eq. © plotted in Fig. [T] demon- 
strates that the approximation is rather accurate within 
the range of tt we are interested in. 
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FIG. 1: Comparison of the left-hand (solid curves) and right- 
hand (dashed) sides of Eq. ||B). The top and bottom panels 
are for charm and beauty respectively. 



Thus, we present the weight factor W 5c Eq. ([5]) in the 
form, 



W £c (r T ) 



1 



2 i 2 

-r T / ri 



(7) 



where r\ = a 2 J(l + a 2 c /2(r 2 J/x ,)), and r\ = b 2 J(l + 

In small-rT approximation the dipole-nucleon cross 
sections in ((3]) have the simple form, CTcc(f*T, £2) = 
C{x2)r\, and (Jccg( r T, £2) — j§C(x2) r T- So, the in- 
tegration can be performed analytically, and we arrive 
at 



S p A{b) = 



where for i = 1 , 2 



T_(6,z)+T+(6,z) 



(8) 



(9) 



Notice that due to color transparency the nuclear 
medium is more transparent than in the Glauber model. 
Moreover, the amplitude Eq. ([9]) does not decrease with 
nuclear thickness exponentially, but as a power. 

Finally, we integrate the attenuation factor Eq. (|SJ) 
squared over the coordinates of the production point and 
arrive at the nuclear ratio, which has the form, 



RpA 



in - ri, 
rir 2 
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where for i = 1 , 2 



1 + ^CT A (b) 



7r 2 

l + ^CT A (b) 
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and 
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With this equations we calculated the nuclear ratio 
Ra/p{v)i Eq. (|TU|) . We rely on the realistic Woods-Saxon 
parametrization for the nuclear density [9J. The results 
at -/s = 200 GeV are depicted as function of rapidity in 
Fig. [U by dotted curve. We see that the steep rise of the 
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FIG. 2: (Color online) Dotted curve presents nuclear suppres- 
sion of J 1^1 as function of rapidity in pA collisions calculated 
with the analytic formula Eq. (|12[) . Dashed curve presents 
more accurate numerical calculations with the full dipole cross 
section Eq. ((2]). Solid curve is corrected for gluon shadowing. 
Data are for d-Au collisions at y/s = 200 GeV [lOJ . 

break-up cross section cr Ec (rr, E Ec ) with energy (it triples 
from y — to y — 2) explains well the observed rapid- 
ity dependence of nuclear suppression. The calculations 
should not be continued far to negative rapidities, since 
the regime of long coherence length breaks down there. 
Besides, additional mechanisms, which cause a nuclear 
enhancement at negative rapidities, must be added. 

We also tested how accurate is the result of analytic 
integration Eq. (|10[) and performed numerical calcula- 
tion using the full equation ([3]) with the dipole cross sec- 
tion Eq. @. The results is plotted by dashed curve in 
Fig. [21 It is pretty close to the previously calculated dot- 
ted curve. Both agree with data within large errors and 
normalization uncertainty. 

In Fig. [3]we present the impact parameter dependence 
of nuclear suppression for J produced with different 
rapidities in p-Au collisions at RHIC. As expected, the 
strongest dependence on rapidity comes from most cen- 
tral collisions. 

We also performed calculations for nuclear effects in 
T production. The only difference with charmonium 
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FIG. 3: 6-dependence of the nuclear ratios for J/\P produced 
with rapidities y — 0, 1, 2, 3 in pAu collisions at y/s = 
200 GeV. 



production is the heavier quark mass, m\, = 4.5 GeV 
and new values of parameters in the parametrization 
Eq. dU, ab = 0.162, bb = 0.037. How accurate is this 
parametrization for bottom quarks is demonstrated in 
the bottom panel of Fig. [TJ The results for nuclear sup- 
pression of T produced on lead at the energies of RHIC, 
y/s = 200 GeV, and LHC, y/s = 5.5 TeV, as function 
of rapidity are depicted in Fig. 21 The upper solid and 

1 1 i 




FIG. 4: Nuclear suppression of T production as function of 
rapidity in pA collisions at y/s = 200 GeV (upper curve, gold) 
and y/s = 5.5 TeV (two bottom curves, lead). The bottom 
dashed curve includes only the effects of bb dipole break-up 
and high twist shadowing of beauty production, solid curve is 
corrected for gluon shadowing. 

bottom dashed curves are calculated at y/s — 200 GeV 
and 5.5 TeV respectively. Notice that in the former case 
the "frozen" approximation is not well justified at y = 0, 
where t c ~ 4 fm. 



B. Leading twist gluon shadowing 

In terms of the Fock decomposition gluon shadowing 
for cc production is associated with higher Fock states 
\ccg) , etc. First, one should evaluate the kinematic con- 
dition for gluon shadowing, t c / 9 > Ra, where t^ cg is the 
coherence time, or the lifetime of a ccg fluctuation in a 
gluon. This time can be related to the Ioffe time as, 

tf 9 = (13) 
xrriN 

where the factor P g « 0.1 was evaluated in [ll[ and found 
to be scale-independent. Its smallness is caused by the 
large intrinsic transverse momenta of gluons in hadrons, 
supported by numerous evidences in data [IH, [l3[. Thus, 
shadowing for gluons onsets at x < 0.01, which is a small 
x value than for quarks. As the result, no gluon shadow- 
ing is possible for charmonium production in any of fixed 
target experiments performed so far 0, [l4[ • 

Even at the energy y/s — 200 GeV the values of X2 
defined in Eq. ^ are too large for gluon shadowing, 
0.024 > X2 > 0.0033 within the measured rapidity in- 
terval < y < 2. We rely upon the NLO analysis [15[ of 
DIS data, which suggests a very weak gluon shadowing, 
as is depicted in Fig. (6) of [l5|. Such a weak shadowing is 
in good agreement with the theoretical predictions |12l |. 
The nuclear ratio presented by dashed curve in Fi g. [2] 
corrected for gluo ns shadowing at Q 2 = 10 GeV 2 [15| 
(see Fig. 6 in [ljf), is depicted by solid curve. We see 
that the effect of gluon shadowing is indeed vanishingly 
small. Even at the energy of LHC, y/s = 5.5 TeV and 
y = gluon shadowing according to [H, Ell is extremely 
small, only 3% at x 2 = 5.5 x 10~ 3 , and will be neglected 
in what follows. 

Apparently, for T production gluon shadowing is 
weaker than for charmonia. No gluon shadowing affects 
T production at the energies of RHIC, as is depicted by 
the upper curve in Fig. [4] since the values of x 2 are too 
large. The effect of gluon shadowing at the energy of 
LHC is is visible and the full calculation is presented by 
the bottom solid curve. It was calculated using the gluon 
nPDF of [H at Q 2 = 100 GeV 2 . 



III. TRANSITION FROM pA TO AA 

At fist glance one might extrapolate the nuclear ef- 
fects from pA to AA collisions in a straightforward way: 
RAA(b, f) = R p a(t) x R p A(b — t), where b is the impact 
parameter of nuclear collisions, and J/^P is produced at 
impact parameter f. Indeed, such a "data driven" pro- 
cedure was used in [Io|, HH to predict the cold nuclear 
matter effects in nuclear collisions based on the mea- 
surements of 6-dependence of nuclear suppression in pA. 
There are several reason, however, which make such a 
transition model dependent. 
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A. Double-color-filtering 

As was discussed above, the survival probability of a 
cc dipole propagating through a nucleus is subject to 
color transparency. This effect in AA collisions turns 
out to be non-factorizable. This can be illustrated 
on the following example. Let a cc dipole of trans- 
verse separation tt propagate through a slice of nu- 
clear medium of thickness Ta with survival probability 
S p a(tt) — exp(— Ct\Ta)- To clarify the appearance 
of the effect we rely on the simplified size distribution 
function W{tt) oc exp[— r|,/(r|i)] adjusted to the distri 
bution Eq. ([5]) with (rf, 
attenuation factor takes the form, 



= 0.045 fm 2 . Then the nuclear 



S pA = / d 2 r T W(r T )S pA (r T ) = 



1 



1 + C{r 2 T ) T A 



(14) 



Naively, one could expect Saa (b) = S 2 A . Such a fac- 
torization is valid only in the dipole representation for 
a given dipole size, Saa^t) = S 2 a {tt) = exp(— 2 x 
Ct^Ta)- Factor 2 is introduced because the dipole is 
simultaneously attenuated by both nuclei. Now we can 
repeat the above averaging over dipole size and compare 
the result (left) with the conventional recipe (right), 



Saa 



I 



1 



l + 2C{r* T )T A [1 + C(4) 



(15) 



We see that the two absorption factors are quite different, 
especially for C{t\)Ta > 1. The source of the difference 
is color filtering. Namely, the mean transverse size of 
a cc wave packet propagating through a nucleus is get- 
ting smaller, since large-size dipoles are filtered out (ab- 
sorbed) with a larger probability pj], [HI ■ Such a dipole 
with a reduced mean size penetrates more easily through 
the second colliding nucleus, compared with a pA colli- 
sion. The mutual color filtering makes both nuclei more 
transparent. 

Now we are in a position to perform realistic calcula- 
tions for the nuclear suppression factor in AB collisions. 
Provided that the cc production occurs in the long co- 
herence length regime for both nuclei, the nuclear sup- 
pression factor at impact parameter b reads, 



Rab(Jo) 



1 



Tab (b) 



d 2 i 



T A (T)T B ((b - f) 
(A+-A A )(A+-A B ) 



where 



x In 



^A = 



(1 + A7 + A+)(1 + A+ + A B ) 



A 



(1 + A++A+)(1 + A7+Ag) 



VM C {Ei)T A {r); 
^L C {Eg)T B (b-r)- 



(16) 

(17) 
(18) 
(19) 



A- 



^C(Eg)T B (b-r) 



(20) 



A B 

and E Sc ' are the energies of the cc in the rest frames of 
the nuclei A and B respectively. The result of Eq. (fTfjj) 
is plotted as the upper solid curve in Fig. [5] For com- 




b (fm) 



FIG. 5: (Color online) Effects of double-color-filtering. J/\& 
suppression by ISI effects in Au-Au collisions at ^fs = 
200 GeV as function of b. The upper and bottom pairs of 
curves (solid and dashed) correspond to y = and energies 
\fs = 200 GeV and 5.5 TeV respectively. Solid and dashed 
curves present the results at y = including and excluding 
the effect of double-color-filtering, respectively. The dotted 
curve demonstrates rapidity dependence of the ISI effects at 
RHIC. It is calculated at y = 2 and is to be compared with 
the upper solid curve at y = 0. 

parison, the result of conventional calculations assuming 
simple multiplication of the suppression factors in the 
two nuclei, is depicted by the dashed curve. We see that 
the mutual color filtering makes the nuclei considerably 
more transparent. This effect should be more prominent 
for production of 'F' and \- 

With Eq. (fTo) we can trace the y-dependence of Raa- 
It turns out to be rather weak at the energy of RHIC, 
due to the approximate linearity of y-dependence in pA 
depicted in Fig. [2 which leads to a compensation of the 
nuclear effects in the colliding nuclei. However, at suffi- 
ciently large ?/, say y — 2, the condition of long coher- 
ence length breaks down in one of the nuclei. Then the 
cc dipole size is not frozen by Lorentz time delation, and 
the filtering in this particular nucleus is not effective any 
more. In this case the conventional multiplicative pro- 
cedure is applicable, but the suppression factor in one 
nucleus (high E Sc ) should be calculated differently, for a 
short t c regime. The result of such calculation is plotted 
by the bottom solid curve in Fig. [3] We see that the 
nuclear suppression at y — 2 is stronger than at y = 0. 
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This happens due to disappearance of the double-color- 
filtering effect. 



B. Boosting the saturation scale in A A collisions 

Another mechanism which violates the conventional 
multiplicative procedure for the transition from pA to 
AA collisions, is the mutual boosting of the saturation 
scale in the colliding nuclei . It is controlled by the fol- 
lowing reciprocity equations [l9| . 



where 



)Ib{xb) = %-ot s {Ql A 



QI)xb9n{x b ,Q 2 sA 



Ql)T B ; 

(21) 



37T 2 

Q 2 sa(xa) = —a s {Q 2 sB 



Ql)x A g N (x A ,< 



(22) 

where we consider a collision of two raws of nucleons Ta 
and Tb , and production of a heavy quark pair with frac- 
tional momenta x A and x b relative to the colliding nucle- 
ons. In what follows we consider heavy quarkonium pro- 
duction at forward rapidities relative to the momentum 
direction of nucleus A, which we call the beam. Corre- 
spondingly, the target nucleus is B. The values of xa(b) 
in (|21[) and (|22|) for charmonium production with positive 
rapidity y are calculated as, 



%A(B) 



(23) 



The mean transverse momentum squared of J/ty pro- 
duced in pp collisions at the energies of RHIC and LHC 
are (py,) « 4 GeV 2 and 7 GeV 2 respectively. 

The gluon distribution function gN(x, Q 2 ) contains the 
parameter Qq needed to regularize the infra-red behavior 
adjusting the saturation momentum in pA collision to the 
known value 12011 . 



QIa&Eqq) = V 2 rT <Tdip(rT,EQ Q ) 

= 2C{E QQ )T A (b), 



r T =0 



T A (b) 



(24) 



where parameter C(Eqq) introduced in Eq. (j3|) is well 
fixed by HERA data. Solution of Eq. flM|) leads to 
the infra-red cutoff parameter Qq ss 1.7 GeV 2 , which is 
nearly independent of energy. 

Solution of equations fl2"Tl) - ([2"2")) shows that the modi- 
fied saturation scales Q 2 A (b) ^ n ^B collisions consider- 
ably exceed the conventional scales Eq. (f24f relevan for 
pA collisions. The essential point of this consideration 
that the boosted saturation scale leads to an increase of 
the break-up cross section of a dipole, because the satu- 
ration momentum is directly related to the factor C(E) 
in the dipole cross section, Eq. (|2"4")l . Thus, the break-up 
cross section UdipifT) for a dipole propagating through 
the nucleus B modifies as 



crdip(r T ) v d i P ( r T) = K A a di p(r T ), 



(25) 



K A = 



QIa 

Q 2 



(26) 



sA 



This boosting factor for the dipole absorption cross sec- 
tion in the nucleus B is controlled by the boosted satu- 
ration scale of the nucleus A, which implicitly depends 
on T A , Tb, x Al and xb in accordance with the equations 
(f2Tjl - (f22|) . Correspondingly the boosted break-up cross 
section of the cc dipole propagating through the nucleus 
A is given by ^ ip (r T ) = K B o- dip (r T )- 

We calculated the factors K A and Kb solving 
the reciprocity equations (|2"T1) - (|2"2"|) supplied with the 
MSTW2008 code |H[ for the gluon distributions. The 
results at -y/s = 5.5 TeV are plotted in Fig. [5] as function 
of T A = Tb (central AA collision) for different rapidities 
of the produced charmonium. Solid and dashed curves 




T A =T B (fm- 2 ) 



FIG. 6: (Color online) The boosting factors Ka(xa) (solid 
curves) and Kb(%b) (dashed curves) denned in (|26[) and cal- 
culated with the reciprocity equations (|2ip - (|22| l for i/s = 
5.5 TeV with Xa(b) defined in (|23p . Each pair of curves is 
marked by the rapidity for which it is calculated. 

present K A and Kb respectively 

We see that at the boosting factors are maximal at 
y = 0, where K A = Kb- At forward rapidities x A rises, 
while xb decreases. As a result, the boosting factors 
K A and Kb remain similar, and both are falling with 
rapidity. Eventually, at large rapidity y — 6 the projectile 
x A becomes so large that the shift of the scale does not 
lead to an increase of the gluons density in A, but makes 
it smaller [22j • Therefore the effect of boosting turns into 
a suppression, i.e. the nucleus A in this case becomes 
more transparent, rather than opaque, for J/^P. 

Similar calculations at y/s = 200 GeV depicted in 
Fig. [7] demonstrate a different behavior. Surprisingly, 
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FIG. 7: (Color online) The same as in Fig. [6] but at -y/i 
200 GeV. 



while Kb is falling with rapidity, as expected, Ka is ris- 
ing. This rise of the boosting factor is related to steep 
drop at forward rapidities of the saturation momentum 
in the projectile nucleus A, the denominator of (I26|) . In 
fact, the very existence of a saturation regime at such 
large xq is questionable. 

We see that the boosting effect significantly increases 
the saturation scales in colliding heavy nuclei, especially 
at the mid rapiditiy and high energies. The rapidity de- 
pendence of the boosting factor is more complicated, and 
is related to different variation with scale of the gluon 
distribution function at different values of Bjorken x. It 
worth reminding that what is plotted in Figs. |5] - are 
the ratios. The absolute value of the saturation momen- 
tum in the nucleus B (A) is steeply rising (falling) with 
rapidity. 

As far as the transparency of nuclear medium for heavy 
quark dipoles in the case of nuclear collisions is different 
from one measured in pA collision the simplified multi- 
plicative prescription of for pA to AA transition 
may be quite incorrect. 

In Fig. [5] we demonstrate the strength of the boosting 
effect on the production rate of J/^P in central (b = 0) 
Au-Au collision as function of impact parameter r. The 
dashed curves include the double-color-filtcring effect, 
but exclude the saturation scale boosting, which is added 
to produce the solid curves. The solid curves include the 
effect of boosted saturation scale. The upper and bottom 
pairs of curves correspond to the energies of RHIC and 
LHC respectively. We conclude that J/'F should be sig- 
nificantly stronger suppressed in AA collisions, than usu- 
ally expected extrapolating from pA, and one should not 
misinterpret this suppression as anomalous effect related 
to FSI with the dense medium. Notice that a stronger 
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FIG. 8: (Color online) Effect of the boosted saturation scale 
on the nuclear ratio for J/ty production in central (b — 0) Au- 
Au collisions at y = as function of impact parameter r. The 
upper and bottom dashed curves correspond to ^/s = 200 GeV 
and 5.5 TeV respectively. They are calculated in the same 
way as the solid curves in Fig. [5] The solid curves here are 
calculated with the boosted saturation scale, which makes the 
nuclei more opaque for heavy dipoles. 



suppression of J/ 1 ^ in the ISI stage of the collision com- 
pared to the one used in [2| should result in even smaller 
transport coefficient of the dense medium extracted from 
RHIC data on J/^ production. 



C. Boosted gluon shadowing 

As we have already mentioned, in terms of the Fock 
decomposition gluon shadowing corresponds to multiple 
interactions of higher Fock states, containing gluons. One 
may wonder if the double color filtering effect can affect 
the amount of gluon shadowing in nuclear collisions. The 
answer is not. The lifetime of such a fluctuation produced 
by a nucleon in the nucleus A may be sufficiently long 
only relative to the nucleus B, but is very short relative 
to the parent nucleus A. Therefore no gluonic fluctua- 
tions undergo double color filtering. In terms of Gribov 
inelastic shadowing this means that the diffractive exci- 
tation of the nucleons of A propagate through B inde- 
pendently of the excitations of B propagating through 
A. 

Thus, the gluon shadowing factors factorize in AB col- 
lisions, 



<>,r) = <(r)<(6-f) 



(27) 



Even if gluon shadowing for heavy quarkonium produc- 
tion in pA collision is known, say, from analyses of DIS 
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data |15[, it is integrated over impact parameter, while 
in (|27[) one needs to know its impact parameter depen- 
dence in order to predict gluon shadowing in heavy ion 
collisions. Already this problem is a serious obstacle for 
extrapolation from pA to AA. 

Aside from the ad hoc parametrizations of the b- 
dependence existing in the literature, one has to rely on 
a fully developed theoretical model for gluon shadow- 
ing to predict its 6-dependence. Unfortunately, no satis- 
factory theoretical description, which would work at all 
kinematic regimes, has been developed so far. The most 
rigorous quantum-mechanical treatment of gluon shad- 
owing within the path-integral technique [jjl, HH is the 
lowest order calculation, which might be a reasonable 
approximation only for light nuclei, or for the onset of 
shadowing. Contribution of higher Fock components is 
still a challenge. This problem has been solved only in 
the limit of long coherence lengths for all radiated gluons, 
in the form known as the Balitsky-Kovchegov equation 
(BK) Hi], HE]- Numerical solution of this equation is quite 
complicated and includes lot of modeling [26] . A simpler 
equation, which only employs a modeled shape of the sat- 
urated gluon distribution, was derived in |22j. It leads 
to a gluon distribution in nuclei, which satisfies the uni- 
tarity bound [13]; an d is quite similar to the numerical 
solutions of the BK equation. The equation reads (22| , 



R 



1 



R 2 n 2 n eff 



(l + R g n ) 2 (l + n eff ) 



(28) 



where 



n (E £ c,b) 
n eff (E 5c ,b) 



9C{E EC ) 



T A (b); 



-C(E- cc )r 2 T A (b). 



(29) 



The energy dependent factor C(E Sc ) was introduced 
above in Q. The mean size of a gluonic dipole, 
0.3 fm, is dictated by data [H, [l3j]. We rely on the satu- 
rated shape of the dipole-nucleon cross section with the 
saturation scale Q,v£EpJ = 0.19 GeV x (E 5c /lGeV)°- u , 
fitted to DIS data [12|, [22| ■ One can switch from energy 
to Bjorken x dependence in these equations, using the 
relation Eq. (|23p and replacing s => 2E 5c mN- 

Notice that Eq. (|28l does not contain any hard scale, 
but only a semi- hard one controlled by tq [13J ■ Therefore 
its solution should be treated as the starting gluon distri- 
bution at the semi-hard scale Q 2 « 4/r 2 w Q 2 , where Qo 
was introduced above in (|2T]) - ([22|) . Shadowing for heavy 
quarkonium production should be DGLAP evolved up to 
an appropriate hard scale. 

The effect of boosted saturation scale leads to a modi- 
fication of the factor C(E 5c ) which is different for nuclei 
A and B, 



C(E 5c ) C A ( B )(Ecc) 



K 



MB) 



C(E 5 



(30) 



Solving equation (|28l) with such a boosted saturation 
scale one arrives at a modified nuclear ratio for gluons 



in nuclei A and B, i?^ S ^(i?c C , b), for which the factor- 
ized relation (|2T|) can be used, 



Rf B (b,r) = Rf(r)R^(b-r). 



(31) 



A numerical example for the boosting effect on gluon 
shadowing is depicted in Fig. [9] The dashed and solid 




T A (fm" 2 ) 



ness Ta calculated with equation ([28} at the semihard scale 
Qo. The dashed curve presents gluon shadowing correspond- 
ing to hadron-nucleus collisions. The solid curve include the 
boosting effects specific for central nucleus-nucleus collisions 
(Ta ~ Tb)- Shadowing is calculated at the starting scale Qo 
and should be evolved up to a higher scale. 



curves show the solutions of Eq. ([28)) without and with 
inclusion of the boosting effects, calculated for J/'J' pro- 
duced with y = at ^/s — 5500 GeV. We see that 
the boosting effect is considerable. As was mentioned 
above, shadowing is calculated at the starting scale Qo 
and should be evolved up to a proper higher scale. 



IV. SUMMARY AND PROSPECTIVES 

The dipole formalism based in the universal dipole 
cross section well fitted to HERA data, successfully de- 
scribes the propagation of heavy quark dipoles through 
cold nuclear matter. Our parameter free calculation pre- 
sented in Fig. [2]well agree with data from RHIC for J/ty 
production in d-Au collisions. We do all calculations 
within the "frozen" approximation, assuming that the 
dipole size does not fluctuate during propagation through 
the nucleus. This condition is satisfied provided that the 
coherence time is long, t c 3> Ra, what is well justified at 
the energies of RHIC and LHC. 
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We identify the main source of the observed nuclear 
suppression as a combined effect of break-up of the pro- 
duced cc dipole in nuclear medium and charm quark 
shadowing. Although both are the high twist effects, 
their contribution to the observed nuclear suppression 
considerably exceeds the effect of leading twist gluon 
shadowing. The latter is expected to be rather weak, 
basing on either theoretical predictions 12[ , or the NLO 
analysis of DIS data. 

Even if the nuclear effects for heavy quarkonium pro- 
duction are understood, the transition to AA collisions 
is rather complicated. Several new phenomena, specific 
for nuclear collisions, make such a transition model de- 
pendent. The first one, double color filtering, makes the 
nuclei more transparent for quark dipoles than in pA col- 
lisions. The second one is the mutual boosting of the 
saturation scales in the colliding nuclei, which makes the 
nuclei more opaque for quark dipoles. Although these 
two effects act in opposite directions, the latter is much 
stronger, as one can see in Fig. [5] 

Another direct consequence of the increased saturation 
scale in nuclear collisions is boosting of gluon shadowing, 
which violates the factorized relation, Eq. (|2?|) . Gluon 
shadowing 

A direct way to access the saturation scale experimen- 
tally is to measure transverse momentum broadening of 
charmonia produced on nuclei. Then the boosting ef- 



fect should show up as an increase of broadening of J/\E' 
produced in AA compared to pA collisions. Indeed such 
an effect was clearly observed in the NA60 and NA50 
experiments at — 158 GeV (28|. The broadening 

in nuclear collisions was found twice as large as in pA 
measurements aith the same path length in the nuclear 
medium. Within the rather good statistical and system- 
atic accuracy, the effect is quite certain. The magnitude 
of the observed boosting is larger than is predicted by 
the equations ([2T]) - (|22[) at this energy, so this problem 
needs further study, and the data are to be confirmed by 
independent measurements. 

Although we considered here only the high energy limit 
of "frozen" dipoles, the pA to AA transition is not trivial 
at medium high energies, say at SPS, either [3]. This 
case will be also investigated and published elsewhere. 
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